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ABSTRACT 



We study the color-magnitude red sequence and blue fraction of 72 X-ray selected galaxy clusters at z = 0.04 - 0.07 from the WINGS 
survey, searching for correlations between the characteristics of the red sequence and the environment. We consider the slope and 
scatter of the red sequence, the number ratio of red luminous-to-faint galaxies, the blue fraction and the fractions of ellipticals, SOs 
and spirals that compose the red sequence. None of these quantities correlate with the cluster velocity dispersion. X-ray luminosity, 
number of cluster substructures, BCG prevalence over next brightest galaxies and spatial concentration of ellipticals. Instead, the 
properties of the red sequence depend strongly on local galaxy density. Higher density regions have a lower RS scatter, a higher 
luminous-to-faint ratio, a lower blue fraction, and a lower spiral fraction on the RS. Our results highlight the prominent effect of the 
local density in setting the epoch when galaxies become passive and join the red sequence, as opposed to the mass of the galaxy host 
structure. 

Key words. Surveys - Galaxies : Clusters : General - Galaxies: evolution - Galaxies: star formation - Galaxies: structure 



1. Introduction 

The location of galaxies in a color-absolute magnitude diagram 
has traditionally been one of the chief methods to examine the 
stellar population properties of galaxies and their evolution. 

The importance of the color-magnitude diagram in 
galaxy studies can ha r dly be overstated. Very early studies 
(Ide VaucouleursI Il96lt IVisvanathan & Sandagd Il977h recog- 
nized the existence in such a diagram of a "red sequence" (here- 
after, RS), where early-type galaxies lie, also known as the 
"color-magnitude relation" (CMR) because it represents a posi- 
tive correlation between galaxy color and luminosity. 

[ The RS represents an easily recognizable feature in the color- 
magnitude diagram that can be used, to first approximation, to 
separate red, passively evolving galaxies devoid of star forma- 
tion, from bluer star-forming galaxies. 

Innumerable studies have investigated the origin of such a re- 
lation, and the reasons why a galaxy belongs to it, or falls below 
it at bluer colors. 

The slope, scatter and location of the RS have long been used 
to place cons traints on the formation epoch of stars in early- 
type g alaxies dBower et al ]|l992l:lEms et al.ll997tlKodama et al.l 
Il998h . The RS is most conspicuous in galaxy clusters, that 
are rich in passive early-type galaxies. A well-d efined RS has 
been observed in clusters up t o high redsh ifts (iLidman et al.l 
[2008; Mei et al. 2009; Strazzu llo et al.ll2010h and the fact that 
the galaxies with the oldest stars preferentially inhabit clusters at 
any epoch can be used to identify galaxy clusters, using their RS, 



up to high redshifts ("Gladders & Yeell 2005HMuzzin et al.ll2009t 
Wilson et al. 2009; Gilbank et al. 20lE. 

Observations of the relative fractions of galaxies on and be- 
low the RS have unc overed evolution in the star format ion ac- 
tivity within clusters (iButcher & Oemleij ll984; Ellings on et al 



'2001 
i2004l: 



iLoh et al ]l2008h an d, more recently, in the field dBell et al 



Faber et all 120071; iBell et al.ll2007l) . It has thus become 



clear that a large number of galaxies have stopped forming stars 
and have turned from blue to red at z below 1 in all environ- 
ments, but in a way that depends on environment. It is also now 
well established that this progressive "passivization" process of 
galaxies proceeds in a downsizing fashion, with more massive 
and luminous galaxies reac hing the RS at ear lier epochs than 
lower-mass, fainter galaxies (ICowie et al.lll996h . Large spectro- 
scopic surveys have confirmed and placed on a very solid sta- 
tistical ground earlier results on the correspondence between 
the bimodality in colors (red and blue) and other main galaxy 
characteristics such as morphological type (early and late) and 
galaxy st ellar mass ( IStrateva et al.ll200"lt iKaufimann et al.ll2003l; 
ICheng et al. 2011 ). 

One of the main outstanding questions remains the phys- 
ical reason why galaxies stop forming stars and become red. 
Naturally, there can be multiple physical mechanisms responsi- 
ble for the end of star formation. Their relative importance might 
vary with galaxy mass, galaxy location (environment), and red- 
shift, among other things (Peng et al. 2010). 

Two mechanisms in particular have been the subject of ex- 
tensive theoretical investigation in the last years: AGN quench- 
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ing, and strangulation by which the hot gas reservo ir is removed 
from galaxies upon merging with a larger halo ( Crot on et al 
2006; Bower et al. 2006; van den Bosch et al. 2008; Guoetal 



201 It iMcCarthv et al.l l2010). Models however are currently un- 



able to reproduce the observed t rends, and generally overpro- 
duce the fraction of red galaxies dFont et al.ll2008l; iBalogh et"al] 
l2009HKimm et al.ll2009h . In clusters, additional mechanisms are 
expec ted and observed to take place, such a s ram pressure strip - 
ping (iGunn & GotlllI972l) and harassment dMoore et al.lll998h . 
though they may be much more efficient than originally thought 
even in lower mass systems than clusters (Bekki 2009). 

The main issue is clearly to what extent the end of star forma- 
tion is due to internal galaxy properties, and to what extent it is 
related to the global environment (for example the halo mass of 
the galaxy host structure) or the local environment. In this paper 
we investigate the relation between the RS properties of galaxies 
in clusters and the global and local environmental properties in 
the local Universe. We use a survey of 77 X-ray selected nearby 
galaxy clusters (the Wide-field Nearby Galaxy-cluster Survey, 
WINGS) to search for correlations between the RS characteris- 
tics and the galaxy environment, with the aim of shedding some 
light on the mechanisms that result in the population of passive 
red galaxies in clusters today. WINGS clusters span a wide range 
of cluster masses, and their galaxies are found in widely differ- 
ent local environmental densities, thus allowing a broad investi- 
gation using a homogeneous, high quality photometric dataset. 

In the following, we describe the WINGS dataset (§2) and 
the method used to define the RS (§3). In §4 we investigate if the 
RS characteristics (slope, scatter and value of the best fit line at 
My - -20), the ratio of the number of luminous to faint galax- 
ies on the RS, the blue galaxy fraction and the morphological 
fractions of galaxies on the RS depend on general cluster prop- 
erties (cluster velocity dispersion. X-ray luminosity, BCG domi- 
nance, number of substructures and concentration of ellipticals). 
We then study the same characteristics of the RS as a function 
of local galaxy density (§4.5). Our conclusions are summarized 
in §5. 

Throughout this paper we use the cosmological parameters 
{Ho, Q„„ Q^) = (70 km s"' Mpc"', 0.3, 0.7). 



2. The Data 

The galaxies examined in this p aper are part of the Wide-field 
Nearby Galaxy-clusters Sun'ey ( iFasano et al. I l2006h . WING^l 
is a multiwavelength survey designed to provide a robust char- 
acterization of the photometric and spectroscopic properties of 
galaxies in nearby clusters, and to determine the variations of 
these properties as a function of galaxy mass and environment. 

Clusters were selected in the X-ray from the ROS AT 
Brigh test Cluster Sample and its extension dEbeling et alJI 19981 
12000) and the X- ray Brightest Abell-type Cluster sample 
dEbeling etanil996l) . WINGS clusters cover a wide range of ve- 
locity dispersion ctcIus, typically between 500 and 1100 kms"', 
and X-ray luminosity Lx, typically 0.2 - 5 x 10'*'*erg/s. 

The survey core dataset, consisting of optical B and V 
imaging of 78 ne arby (0.04 < z < 0.07) galaxy-clusters 
dVarela et al.ll2009l) . has been complemented by several ancil- 
lary projects: (i) a spectroscopic follow up of about 6500 galax- 



ies in a subsample of 48 clusters, obtained with the spectro 
graphs WYFFOS@WHT and 2dF@AAT dCava et 



e s pectro - 
aH 120091) : 



(ii) near-infrared (J, K) imaging of a subsample of 28 clus- 
ters obtained with WFCAM@UKIRT (IValentinuzzi et al.ll2009l) 
; (iii) U broad- and H„ narrow-band imaging of subsamples of 
WINGS clusters, obtained with wide-field cameras a t different 
telescopes (INT, LBT, Bok, see lOmizzolo et al.ll2bl ll in prepa- 
ration). 

The B and V photometric catalogs used in this paper are de- 
scribed in detail in ( Varela et a l. 2009). Briefly, our catalogs are 
90% complete at V ~ 21.7, and the star-galaxy separation was 
visually checked, making the number of misclassifications neg- 
ligible, down to y = 22. The distance moduli of our clusters are 
in the range 36.1-37.4, therefore our photometric catalogs are 
highly complete and reliable down to My = -15.7 or fainter 
The photometry was performed on images in which large galax- 
ies and halos of bright stars were removed after modeling them 
with elliptical isophotes. This procedure greatly improved the 
photometry of the large galaxies and incre ased by 16% the d e- 
tection rate of objects projected onto them dVarela et alJl2009l) . 

WINGS galaxy morphologies were derived from V images 
using the purposely devised tool MORPHOT (Fasano et al., 
2011 submitted, and Appendix A in Fasano et al. 2010). Our 
approach is a generalization of t he non-par ametric method p ro- 
posed bv IConselice et al] d2000l) (see also. IConselicell2003h . In 
particular, we have extended the classical CAS (Concentration/ 
Asymmetry/clumpinesS) parameter set by introducing a num- 
ber of additional, suitably devised morphological indicators, us- 
ing a final set of 11 parameters. A control sample of 1,000 vi- 
sually classified galaxies has been used to calibrate the whole 
set of morphological indicators, with the aim of identifying the 
best sub-set among them, as well as of analyzing how they de- 
pend on galaxy size, flattening and S/N ratio. The morphological 
indicators have been combined with two independent methods, 
a Maximum Likelihood analysis and a Neural Network trained 
on the control sample of visually classified galaxies. The final, 
automatic morphological classification combines the results of 
both methods. We have verified that our automatic morphologi- 
cal classification reproduces well the visual classification by two 
of us (AD and GF). In particular, the robustness and reliability 
of the MORPHOT results turn out to be comparable with the 
typical values obtained comparing the visual class ifications ob- 
tained by different experienced human classifiers dFasano et al.l 
1201 iL in preparation). 

All the following analysis is carried out for galaxies within 
7^200 /2, where /?200 (usually considered an approximation for the 
cluster virial radius) is the radius enclosing the sphere with inte- 
rior mean density 200 times the critical density of the Universe at 
that redshift. The /?20o/2 values were measured from the cluster 
velocity dispersions given in Table 1 as 



■^^200 - 1-73 



1 



1000 km s 1 ^Q^ + Qo(l +2): 



: /r' Mpc 



(1) 



' Please refer to WINGS Website for updated details on the survey 
and its products, http : //web . oapd. inaf. it/wings 



After careful visual inspection of all the B-V versus V color- 
magnitude diagrams (CMDs) of our clusters, we decided to in- 
clude only 72 clusters in this work because of incomplete radial 
coverage out to /?20o/2 or nearly absent red sequence (RS). 

In Table[T] we present a number of global properti es of our 
clusters. Cluster redshifts are from ICava et all d2009l) . Cluster 
velocity dispersions were derived, following the recipes given 
in Cava et al. (2009), from the WINGS database that collects, 
together with our data, all the available redshifts from NED. X- 
ray total luminosities are taken from Ebeling et al. (1996, 1998, 
2000), and they have been converted to our adopted cosmology. 
The number of substructures found in each cluster are drawn 
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from the analysis of lRamella et alJ (l2007h . who searched for sub- 
structures in the projected distribution of galaxies in WINGS im- 
ages using an adaptive-kernel procedure. The BCG prevalence 
value is determined with the following equation: 



BCGpreV = (Vrank2 + Vrank3)/2 - V^ankl 



(2) 



where the V magnitudes refer to the first (the BCG), second and 
third ranked galaxies of the cluster The higher the BCG preva- 
lence value, the greater is the separation in magnitudes between 
the BCG and the next brightest cluster galaxies. It is believed 
that a system of galaxies, where most of the mass has been as- 
sembled very early, develops a larger magnitude gap between 
the brightest and second brightest galaxies compared to sys- 
tems that form later (Dariush et al. 2010 and references therein), 
therefore the BCGprev value would be somehow related to the 
main cluster "assembly epoch". The concentration of elliptical 
galaxies (EsConc) is measured through a modified versio n of 
the Gini coeflicient dXbraham et al. 2003i iLotz et al.ll2004 . In 
particular, we define this coefficient as the diff'erence between 
the area subtended by the cumulative distribution function of 
the elliptical galaxies (rank-ordered according to their cluster- 
centric distance) and the diagonal of the s quare (Lorentz curve ; 
for a perceptual i llustration see figures 1 in lAbraham et al.l2003t 
ILotz et all l200l . In this formulation, the Gini coefficient in- 
creases as the concentration of galaxies toward the cluster center 
increases. 

Finally, we measured galaxy local densities for each galaxy 
in our sample using the circular area (Aio) containing the 10 
nearest projected neighbors in the photometric catalog (with 
or without spectroscopic membership) whose V-band absolute 
magnitude would be < -19.5 if they were cluster members. As 
we only want to count as neighbours the members of the clus- 
ter, a statistical field correc t ion ha s been applied to the counts 
using Table 5 in iBerta et al.l ( |2006|) . In particular, since the field 
counts in the area containing the 10 nearest neighbors are not 
integer numbers, Aio is obtained interpolating the two A„ ar- 
eas for which the corrected counts (or the number of spectro- 
scopic members, if greater than them) are immediately lower and 
greater than 10. A similar interpolation technique has also been 
used when the circular area containing the 10 nearest neighbors 
is not fully covered by the available data (galaxies at the edges 
of the WINGS field). In this case, at increasing n (and the cor- 
responding area A„), a coverage factor has been evaluated as the 
ratio between the circular area and the area actually covered by 
the observations. Then, the counts n have been scaled upwards 
to account for the corresponding coverage factors and then have 
been corrected for the field counts. Finally, as in the previous 
case, Aio has been obtained by interpolating the two A„ areas 
for which the corrected counts are immediately lower and greater 
than 10. 

3. Red Sequence and Field Subtraction 

One of the main issues in defining the red sequence of our clus- 
ters is the removal of the interloper galaxies. For all galaxies 
with redshift information, either from WIN GS spectra or lit era- 
ture data, the cluster membership is known dCava et alj2009l) . In 
the absence of spectroscopy, the only alternative way is a statis- 
tical approach. 

3. 1 . The interloper issue 

The Monte Carlo statistical field subtraction tecnique is 
the standard method used (see e.g. iKodama & Bowe 



iPimbblet etani2002h . We used the iBerta et al.l (l2006h catalog 
of field galaxies in the ELAIS-Sl area (hereafter simply called 
FIELD), a cluster-free portion of the sky, to determine the num- 
ber of contaminating galaxies as a function of magnitude and 
color Both the WINGS (which comprehend cluster+field galax- 
ies) and the FIELD data have been conveniently binned onto a 
grid in the color-magnitude diagram space. Obviously, for each 
cluster, the numbers of field galaxies in the grid are rescaled to 
the size of the field covered by WIN GS. 

Following lPimbblet et alJ(l2002h . we defined the probability 
for each galaxy to be an interloper as follows; 



Pinticol, mag) = 



Np{col, mag) 
N-wicol, mag) 



(3) 



where N-p and A^w are the number of FIELD and WINGS galax- 
ies in each (col, mag) bin, respectively. We further expanded the 
method to include all the redshift information we have: spec- 
troscopically confirmed members A^yes and non-members A^no 
were always retained and rejected, respectively, in the analysis. 
Eq.© was modified accordingly: 



Np(col, mag) — N^oicol, mag) 



Nwicol, mag) - Nno{col,mag) - Nyes{coI, mag) 



(4) 



One of the main caveats of this method is the possibility that 
the proba bility becomes negativ e: in this case we used the same 
recipe of iPimbblet et al.l (l2002h by applying an ad aptive grid 
tecnique (for further details see IPimbblet et al.ll2002l Appendix 
A). Once the interloper probability was determined at each grid 
position, we performed 100 Monte Carlo simulations generating 
a number between 0.0 and 1.0 for each galaxy. If the random 
number was greater than the corresponding Pint, than the galaxy 
was flagged as a member, otherwise it was assigned to the inter- 
loper population. 

With this procedure the assumption that the FIELD sample 
is a perfect realization of the real interloper population is miti- 
gated. Indeed, it can be demonstrated that with this procedure we 
are actually giving to each bin of the FIELD with galaxies, a 
corresponding poissonian error of ; thus, it is like extracting 
at each run a new FIELD sample (a sort of bootstrap). 

In FiglT] we show the 37-th run of a Monte Carlo realiza- 
tion of the cluster Abell 119. In the upper panel all galaxies 
from the catalog are presented, but shown as green circles and 
blue crosses are the spectroscopic members and non-members, 
respectively. The use of the spectroscopic membership informa- 
ti on is the only d i fferen ce from the interloper subtraction used 
m IPimbblet et all (|2002|) . Red open circles and black filled dots 
refer to the galaxies assigned to the non-member and member 
population by the simulation. We also show the different mor- 
phologies, when available (see legend in the figure), of all cluster 
members (central panel) and only RS galaxies, as defined in the 
next section (bottom panel). Note the non-negligible presence of 
red spirals (blue and cyan open circles). 

3.2. Defining the red sequence 

We calculate rest frame absolute magnitudes and colors by ap- 
plying a galaxy-type dependent k-correction based on the ob- 
served (B-V) color of the galaxies within an aperture (diam eter) 
of lO.Skpc using the K-corrections from (iPoggiantill 19971) . The 
physical aperture assures a well-sampled, distance-independent 
color for all galaxies. We further apply a cut in color 0.5 < 
(B - V)ioMpc ^ 1.5 to eliminate those eventual (few) galaxies 
that are surely either too red or too blue to belong to our clusters 
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Fig. 1. Visual example of a Monte Carlo realization for the cluster Abell 1 19. In the upper panel, the red empty dots and the black 
filled dots are the non-members and members assigned by the simulation (run 37th), respectively. Shown as green circles and 
blue crosses are the spectroscopic members and non-members, respectively. In the middle panel only members according to the 
Montecarlo realization are shown, but divided in different morphological types (see legend and the figure in electronic format). In 
the bottom panel, the RS extracted is shown, again with the morphological color coding. The dotted vertical lines mark the absolute 
V magnitude limits used to define the Luminous and Faint population on the RS. 



but might have remained after the statistical subtraction. In the 
following, we consider only galaxies brighter than My ~ - 1 8 
for which both our photometric and morphological catalogs are 
complete (Varela et al. 2009, Fasano et al. 201 1). 



We use the robust bi-weight estimator fitting procedure in 
two steps to identify the RS. The first step consists of fitting all 
the galaxies from the previous selection process with absolute 
V magnitude -21.5 < My < -18 and calculating the median 
distance dred of all the galaxies from the best fit. We choose this 
magnitude interval because most of our CMDs (see FiglT]i ex- 
hibit a change of slope (are shallower) at higher luminosities. 



of the second fit are assigned to the RS-populationQ All non RS- 
galaxies bluer than that are assigned to the blue-population. 

This procedure is applied to each of the 100 Montecarlo 
Simulations and a final mean value and its error are found for 
each individual cluster (see Tabl^H). In Fig.© we present a sam- 
ple plot of one of the Montecarlo realizations for each cluster, 
with the best fit parameters and the global scatter for the red 
sequence and the main properties of the cluster shown at the 
bottom of each panel. The vertical lines indicate the V-band ab- 
solute magnitude limits used to define the "Luminous" and the 
"Faint" red sequence population (see section |4|2]). 



In the second step only galaxies with -21.5 < My < - 1 8 and 
distance 3 xdred from the previous step are fitted, and the RMS of 
the residuals is calculated. Only galaxies within 0.2 magnitudes 



^ Using galaxies within 2 x RMS from the fit did not change any of 
the results of this paper. 
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Fig. 2. A random realization of the Montecarlo interloper subtraction technique. Only galaxies assigned to the clusters are shown: 
red dots are morphologically early-type galaxies (Es and SOs), blue dots are late-type ones and black dots are galaxies with no 
available morphological classification (colors only available in the on-line version). The best fitted RS of the 100 realizations is 
drawn (solid black line) with ±0.2mag limits (dashed black lines) identifying the red sequence. The RMS numerical values are 
reported at the top of each panel. The global properties of each cluster are reported, for reference, at the bottom of each panel. The 
vertical lines mark the absolute V magnitude limits used to define the Luminous and the Faint population of the RS. 
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Fig. 3. Slope, value of best fit line at My - -20 and scatter of the median RS values of the 100 Montecarlo runs for each cluster 
versus cluster velocity dispersion {(TcIus, left panels) and cluster X-ray luminosity {Lx, right panels). Error bars refer to the error of 
the median of the 100 simulations. In red (on-line version), clusters whose RS is contaminated by another structure (see text). 



4. Results 

4.1. Slope and scatter 

In Tabl^Il in addition to the global parameters of our cluster 
sample, we present the median slope and scatter of the 100 
Monte Carlo realizations described in section 13.21 The median 
slope for all our clusters is -0.047 + 0.001, and the typical scat- 
ter is of the order of ~ O.OSmag. 

In Figl3]we plot the slope, the fit value at My - -20 and 
the scatter of the RS against the cluster velocity dispersion and 
total X-ray luminosity (both of which are related to the cluster 
mass). We have coloured in red those clusters whose RS appears 
to be contaminated by another galaxy structure visible in Fig. 2 
(A133, A151, A548b, A2149, A2382, A2415, A3490, A3667, 
A3809)3 although it should be kept in mind that several other 
clusters have a very broad RS suggestive of more than one struc- 
ture along the line of sight and close in redshift to the main clus- 
ter. 



' We have not highlighted here those clusters with a second RS that 
lies at sufficiently higher redshift not to be expected to influence the RS 
of the main WINGS structure, such as A2572, A3560, A3490, RX1022. 



If the epoch and rate at which galaxies turn red and reach the 
color-magnitude RS depended upon the mass of the cluster to 
which galaxies belong at z = 0, we should observe a correlation 
between cluster mass and RS properties. Instead, we do not find 
any such correlation in our sample of 72 local clusters. 

Similarly, we do not find any dependence of the RS slope, 
scatter, and My = -20 value on the number of substructures, 
the BCG prevalence value and the concentration of ellipticals 
(not show n ). Ou r results agree with the previous findings of 
IStott et aP (l2009h which showed a lack of correlation between 
the RS slope and the cluster cr, Lx and BCG dominance at 
z = 0.1-1. 



4.2. Luminous-to-faint ratio and biue fractions 

The number ratio of luminous-to-faint RS galaxies is often used 
to investigate the evolution of the faint end of the red lumi- 
nosity function of cluster galaxies. High-z studies find a de- 
ficiency of faint red galaxies in comparison to local galaxy 
clusters, as the luminous-to-faint ratio is seen to evolve from 
hig her values at high-z to l ower values at low-z ( see, among oth- 
ers, iDe Lucia et all (l2004 : iTanaka et al] (l2005h : iDe Lucia etaP 
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C clus 
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slope 


IviVio 


A85 


0.052 


1052.0 


44.92 
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1.405 


0.236 


-0.039 ± 0.007 


0.024 


A119 


0.044 


862.0 


44.51 
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0.990 


0.179 


-0.042 ± 0.007 


0.025 


A133 


0.060 


810.0 


44.55 




0.575 


0.067 


-0.030 ± 0.011 


0.043 


A147 


0.045 


666.0 


43.73 


' 


0.505 


0.162 


-0.043 ± 0.012 


0.025 


A151 


0.053 


760.0 


44.00 


2 


1.220 


0.108 


-0.051 ± 0.007 


0.047 


A160 


0.044 


561.0 


43.58 


3 


0.485 


0.119 


-0.037 ± 0.010 


0.044 


A168 


0.045 


503.0 


44.04 


1 


1.130 


0.128 


-0.064 ± 0.013 


0.066 


A193 


0.049 


759.0 


44.19 





1.000 


0.164 


-0.046 ± 0.011 


0.022 


A376 


0.048 


852.0 


44.14 


1 


1.345 


0.165 


-0.041 ± 0.008 


0.031 


A500 


0.068 


658.0 


44.15 
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0.435 


0.175 


-0.037 ± 0.009 


0.033 


A548b 


0.044 


848.0 


43.48 




0.635 


0.127 


-0.042 ± 0.008 


0.033 


A602 


0.062 


720.0 


44.05 


1 


0.595 


0.118 


-0.029 ± 0.011 


0.047 


A671 


0.051 


906.0 


43.95 


2 


1.070 


0.207 


-0.034 + 0.008 


0.038 


A754 


0.055 


1000.0 


44.90 
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1.140 


0.164 


-0.042 ± 0.006 


0.030 


A780 


0.057 


734.0 


44.82 




1.220 


0.111 


-0.043 ± 0.011 


0.030 
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0.045 


710.0 


43.89 




1.680 


0.160 


-0.037 ± 0.013 


0.030 


A970 


0.059 


764.0 


44.18 
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0.290 


0.214 


-0.044 ± 0.010 


0.029 


A 1069 


0.065 


690.0 


43.98 
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1.030 


0.186 


-0.042 ± 0.011 


0.043 


A1291 


0.051 


429.0 


43.64 
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0.144 


-0.044 ± 0.014 


0.027 
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0.046 


640.0 


43.86 





1.385 


0.098 


-0.059 ± 0.009 


0.046 


A 1644 


0.047 


1080.0 


44.55 
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1.275 


0.000 


-0.035 + 0.006 


0.028 


A 1668 


0.063 


649.0 
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1.060 


0.153 


-0.059 ± 0.014 


0.038 
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0.046 


853.0 
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0.145 
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-0.044 ± 0.007 


0.043 
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0.063 


725.0 


45.05 


1 


1.480 


0.185 
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0.032 
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0.063 
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0.046 
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0.045 
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0.028 


A1991 
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0.855 


0.102 


-0.048 ± 0.009 


0.037 


A2107 


0.041 


592.0 


44.04 





1.395 


0.080 


-0.025 ± 0.011 


0.028 


A2124 


0.067 
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44.13 
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0.124 


-0.050 ± 0.008 


0.039 
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0.124 
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IGilbank et all (l2008h: iLu et al 
('2009V. 'C apozzi et al 



2010) and Andreon (2008); Crawford et al. (2009) for opposite 
conclusions). This is generally interpreted as evidence for a large 
number of relatively faint galaxies having moved to the RS only 
recently, most probably due to the quenching of the star forma- 
tion caused by the high-density environment. The luminous-to- 
faint ratio has been found to evolve in the field too, but to be 
lower in clusters than in the field at all redshifts out to z = 1, 
implying that t he faint end of the red seq uence was established 
first in clusters (IGilbank & Baloghll2008l) . 

We wish to study how the luminous-to-faint RS ratio changes 
from cluster to cluster and depends on cluster properties at low 
redshift. To define the two galaxy populations we corrected 
our absolute magnitudes for passive evolution: we consider all 
galaxies brighter than My = -20 - lumcon to be luminous, 
and all galaxies fainter than this magnitude and brighter than 
My = -18.2 - lumcon to he faint, where lumcon is the pas- 
sive evolution correction with respect to z = 0. We recall here 
that these are the same absolute V magnitude limits adopted by 
iDe Lucia et a l. (2007) and most other works at high redshift, but 
that high-z studies usually consider the (U-V) rest-frame color, 
therefore a direct comparison cannot be performed. 

Our large sample of local clusters show a large scatter in the 
luminous-to-faint ratio. In the top panels of Fig|4]and Fig|5]we 
plot the luminous-to-faint ratio as a function of cluster veloc- 
ity dispersion. X-ray luminosity, number of substructures, BCG 
prevalence value and concentration of ellipticals. 

As was the case for the RS parameters, no correlation be- 
tween the luminous-to-faint ratio and the cluster properties is 
found. This agrees with recent results from lCapozzi et all (1201 Ol) 
who found no relationship between luminous-to- faint ratio and 
cluster X-ray luminosity at low redshift, and with lRudnick et all 
(I2009h who failed to detect a dependence of the red galaxy lu- 
minosity function on cluster velocity dispersion in the SDSS. 

Such lack of a correlation suggests that the downsizing trend 
in star formation (the different distribution of times at which 
galaxies turn red and reach the RS depending on their luminos- 
ity) does not depend in a simple way from global cluster proper- 
ties such as the cluster mass, or level of substructure. However, 
the large scatter of luminous-to-faint values in our sample shows 
that there is a great diversity in the way the RS is populated 
as a function of galaxy magnitude. Most likely, this indicates 
widely different quenching histories in faint galaxies from clus- 
ter to cluster, even for clusters of the same mass at z = 0. 

The blue fraction is defined as the ratio between the num- 
ber of blue galaxies and the total number of galaxies in the 
cluster, for My < -18.2. In the last decade the blue frac- 
tion has been observed to have a large scatter of values at any 
given redshift, and to depend on several obvious parameters 
such as the magnitude limit used and the cluster-centric dis- 
tance. However, contrasting conclusions have been reached as to 



richness and presence 


of substructure (iButcher & Oemlerl 


1984; 


Wans & Ulmeil 


1997 


;ISmailetal. 19981 Metevier et al.l 


200C 


Ellinason et al.| 


2001 


i iMaraoniner et al.' '2001'; 'Pimbblet et al. 



2002; Fairlev et al.l 12002 ; Balogh et al. 2004; De Propris et alj 



2004 lBarkhouseetal.1 120091) . All these findings together im 



ply that environmental effects may compete with, and possibly 
mimic, evolutionary trends. Thus, it is of paramount importance 
to try to disentangle the cluster-to-cluster variance and the de- 
pendence on cluster properties in order to properly understand 
and interpret the studies of high-z clusters and to investigate how 
the current galaxy populations in clusters came to be. 



Figs|4] (bottom panels) shows that in WINGS there is no 
dependence of the blue fraction on either cluster velocity 
dispersion or total X-ray luminosity, as f ound previously in 
the local Universe and at higher redshifts dFairlev et al.l 120021 : 
IDe Propris et al.ll2004t lGotoll200l . and in agreement with a flat 
median star-forming (emission-line) fra ction for clusters with 
cr > 500km s~' in the local Universe (iPoggianti et al.l 120061 : 
iPopesso et alj|20()7l) . Furthermore, we find no trend with the 
number of substructures, the BCG prevalence value and the con- 
centration of elliptical galaxies (Fig|5l). A quite large scatter in 
the blue fraction is observed fixing any of these parameters, 
though most clusters have a blue fraction below 20%, with a 
me dian of 0.16±0.0 3. 

IPoggianti et alJ (l2006l) proposed a scenario in which the pop- 
ulation of passive galaxies (those devoid of ongoing star forma- 
tion at the time they are observed) consists of two different com- 
ponents: primordial passive galaxies (the most massive, mostly 
ellipticals), whose stars all formed at z > 2-3, and quenched 
galaxies (on average less massive, mostly SOs), whose star for- 
mation has been truncated due to the dense environment at later 
times. Comparing with simulations, they found that at z = the 
observed fraction of passive galaxies in clusters resembles the 
fraction (in mass and in number of galaxies) that has resided in 
clusters (M^ys > IO^'^Mq) during at least the last 3 Gyr. 

In this picture, the median fraction of star-forming galax- 
ies, and the median blue fraction we observe, do not correlate 
with the cluster cr or L^.This is because in clusters more mas- 
sive than 500km s"' at z = 0, the median fraction of galaxies that 
have spent enough time in a massive environment to have their 
star formation quenched by some (unidentified) environmental 
mechanism does not vary systematically with cluster mass. This 
may be a viable explanation also for the lack of correlation be- 
tween the luminous-to-faint ratio and the global cluster proper- 
ties, as the blue galaxies and the faint red galaxies together form 
the overall faint cluster population subject to quenching. 

4.3. Morphological fractions and red spirals 

When studying the CMD and, more generally, the evolution of 
galaxy clusters, it is important to evaluate the fractions of the dif- 
ferent morphological types residing on the red sequence. Indeed, 
apart from the well known presence of red Es in galaxy clusters, 
it is not completely clear if there is some correlation between the 
frequency of red SOs and late-type galaxies and the global prop- 
erties of the clusters. The presence of SOs and late-type galaxies 
on the RS is often ascribed to the capacity of the hostile cluster 
environment to quench star formation or even strip and/or induce 
a fading of the disk of galaxies of the later types. 

It is apparent from Fig|6]that there is no dependence of the 
RS morphological fractions on either X-ray luminosity and clus- 
ter velocity dispersion. Similarly, no correlations were found 
with the number of substructures, the BCG prevalence value and 
the concentration of ellipticals (plot not shown). In the scenario 
of a double channel for passive galaxies discussed above, this 
is expected given that the median fractions of both primordial 
(ellipticals) and quenched galaxies (SOs) are flat with the cluster 
velocity dispersion (see Fig. 16 in Poggianti et al. (2006)). 

The median morphological fractions of galaxies on the RS 
are 0.35±0.05 for eflipticals, 0.46±0.04 for SOs, 0.17+0.04 for 
early spirals (Sa to Sbc, < Ttype < 4) and 0.02+0.02 for later 
type spirals (4 < T < 8). 

It is interesting to compare these values with those found by 
ISanchez-Blazquez et al.l ( l2009l) in clusters atz = 0.4-0.8, which 
were obtained adopting similar cluster radial limits and galaxy 
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Fig. 4. Luminous-to-faint ratio (top panels) and blue fraction (bottom panels) vs. cluster velocity dispersion and total X-ray lumi- 
nosity in WINGS clusters. Errorbars are the mean error on the 100 Montecarlo realizations. 



magnitude (passively evolved) limits to our study. These authors 
found an increase in the late-type fraction on the red sequence, 
from 25% at z = 0.75 to 44% at z = 0.45, due to the fact that the 
red-sequence becomes more populated at later times with disc 
galaxies whose star formation has been quenched. Interestingly, 
however, the same authors predict a subsequent decrease of the 
red late-type fraction from z - 0.45 to z = 0, due to the combi- 
nation of the strong morphological evolution from spirals to SOs 
and the mild variation in the RS luminous-to-faint ratio observed 
in this redshift interval. 

Indeed, our ~ 20% fraction of late-type galaxies on the RS is 
much lower than the value observed at z = 0.45, and is consistent 
with a scenario in which over the interval z = 1-0.4 the RS 
accretes spirals that have stopped forming stars, a quite large 
fraction of which are transformed into SOs at z < 0.4 (Desai et 
al. 2007, Dressier et al. 1997, Fasano et al. 2000, Postman et al. 
2005). 



4.4. Elliptical to SO ratios and cluster evolution 

Regarding the formation history of galaxy clusters, one may ex- 
pect that clusters with a large presence of red spirals may also 
have a higher fraction of SOs compared to Es on the RS. In other 
words, clusters with a larger number of passive late-type galax- 
ies could have transformed a larger number of them into SOs. 



Figure|7]shows the relation between the RS late-type fraction 
and the RS ratio of Es to SOs. Clusters with a larger fraction of 
red spirals tend to have a population of early-type galaxies on the 
RS dominated by SO morphologies. According to the Spearman 
test the E/SO number ratio on the RS is anticorrelated with the 
late type fraction on the RS with a 95% probability. 



5. The main driver: the local density 

In the previous sections we have shown that the RS param- 
eters, luminous-to-faint ratio, blue fraction and morphological 
fractions on the RS do not depend on cluster global properties 
such as velocity dispersion or X-ray luminosity, nor on elliptical 
galaxy concentration, BCG prevalence or number of substruc- 
tures. 

In Figsl8]and|9]we show that all the galaxy properties con- 
sidered do, however, vary with the local galaxy density. Most 
noticeably, higher density regions have a tighter RS (lower RS 
scatter), and a lower fraction of blue galaxies. The RS generally 
becomes steeper (RS slope increases) with local density, except 
for a high value in the lowest density bin. We note that a depen- 
dence of the location of the CM relation and of the blue fraction 
on local density has been found before in the Shapley superclus- 
ter at z=0.05 by Haines et al. (2006). 
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Fig. 5. Luminous-to-faint ratio (top panels) and blue fraction (bottom panels) vs. number of substructures, BCG prevalence value 
and concentration of elliptical galaxies in WINGS clusters. Errorbars are the mean error on the 100 Montecarlo realizations. 



The high RS slope in the lowest density bin is accompanied 
by slighly higher values of the fraction of ellipticals and of the 
E/SO number ratio compared to denser regions (Fig|9|l. It is hard 
to assess the significance of these single points, and wider im- 
ages will be needed to sample lower density outer regions of 
clusters. 

Interestingly, the RS luminous-to-faint ratio increases mono- 
tonically with local density (Fig. 8), implying that the relative 
proportion of luminous red galaxies and faint red galaxies is 
higher in denser regions. This is the first time such a trend is 
observed. At high-z, the opposite trend has been found, con- 
sistent with the fact that the build-up of the CMR is delayed 
in low er density environments (iTanaka et alJl2005t iTanaka et al.l 
l2007l) . The increase of the luminous-to-faint ratio with density in 
WINGS clusters seems to reflect, instead, the fact that WINGS 
higher density regions host proportionally more high-mass than 
low-mass galaxies (regardless of color) than lower density re- 
gions, as shown in Vulcani et al. (201 1). 

Finally, higher density regions have a lower fraction of spi- 
rals on the RS than lower density regions (Fig|9l): the RS in the 
highest density regions includes only a few spirals (~ 10%), 
while the RS in low density regions is composed by up to 30% 
of spirals. The relative proportion of ellipticals and SOs does not 
change with local density, except for the highest density regions 
where ellipticals dominate. 



It is important to stress that the correlations between RS 
properties and local density shown in Fig|8] are not driven by 
differences between the RS properties of the different morpho- 
logical types in combination with the morphology-density rela- 
tion. In particular, we wished to test whether a larger fraction 
of spirals at low densities might induce the trends observed. We 
found that none of the conclusions from Fig[8]change if we only 
include ellipticals, only SOs, or only ellipticals+SOs. 

6. Discussion and conclusions 

In this work we have analyzed the properties of the (B-V) versus 
V color-magnitude red sequence of 72 X-ray selected clusters at 
z = 0.04 - 0.07 from the WINGS survey. 

We have searched for correlations between the characteris- 
tics of the RS and the properties of both the global and local 
environment, finding the following main results: 

- The location and parameters of the RS, as well as the magni- 
tude and morphological distributions of its galaxies, do not 
depend on global cluster properties. Specifically: 

- neither the slope nor the scatter of the RS, nor 

- the number ratio of luminous-to-faint galaxies on the RS, 
nor 

- the fraction of blue galaxies, nor 
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- the morphological mix on the RS, that is, the fraction of 
galaxies on the RS that are ellipticals, SOs and spirals, 

none of the above vary systematically with cluster velocity 
dispersion nor X-ray luminosity (hence with cluster mass), 
nor with the number of substructures, BCG magnitude sepa- 
ration from the next brightest galaxies nor spatial concentra- 
tion of elliptical galaxies. 

We find a weak tendency for clusters whose RS is "rich" in 
SO galaxies (compared to the population of ellipticals) to be 
also "rich" in red spirals. The populations of SO and spirals 
on the RS are coupled at some level. 

The strongest correlations we find are between the RS and 
the local galaxy density. The scatter and, possibly, the slope 
of the RS vary with local density, as do the luminous-to- 
faint ratio, blue galaxy fraction and morphological mix of 
galaxies on the RS. The latter means that, remarkably, a clear 
morphology-density relation, especially for spirals, is visible 
even restricting the analysis to galaxies on the red sequence. 

However, the correlations between RS properties and local 
density are not due to the morphology-density relation (i.e. a 
higher fraction of spirals at lower densities), as all the trends 



in Fig|8] persist even when restricting the analysis to only 
ellipticals, only SOs, or only early-type galaxies. 

These trends suggest that galaxies in higher density regions 
within clusters became passive and evolved morphologically to- 
wards earlier types at earlier epochs than galaxies in lower den- 
sity regions that are today within the same clusters, as witnessed 
by the lower RS scatter, lower blue fraction, and lower spiral 
fraction on the RS at higher densities. 

Local density thus appears to be the main factor govern- 
ing the evolutionary pace of galaxies in these clusters. If this 
is the case, the lack of a trend of the RS characteristics with 
the cluster mass is not surprising, given that the distribution 
of local galaxy dens ities does not depend on cluster mass 
( lPoggiantietal.1120101) . 

In the same way as the morp hological mix o f galaxies (of all 
colors) d epends on local density (lDresslerll980h and not on clus- 
ter mass dPoggianti et alj|2b09l) in the local Universe, this work 
highlights the prominent effect of the local density as opposed to 
the mass of the galaxy host structure in setting the passivization 
epoch of galaxies and their arrival on the RS, and the subsequent 
morphological transformation to SOs. 
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The dependence of the RS properties on local density can 
originate from several mechanisms, and at different cosmolog- 
ical epochs. Since the present local density is expected to cor- 
relate with the initial density in the primordial phases, the star 
formation process may be accelerated in dense regions from the 
start, yielding shorter formation times. Physical processes act- 
ing at later times, among those commonly refereed to as envi- 
ronmental mechanisms, can accelerate and/or quench star for- 
mation more efficiently in denser regions. Among these, there 
are gravitational interactions between galaxies, before or after 
the accretion of galaxies onto clusters, and gas removal mech- 
anisms. While a detailed modeling of such processes is beyond 
the scope of this paper, the combination of the lack of trends with 
cluster mass and of the observed correlations with local density 
represent a solid observational contraint for future studies. 
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